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^ '. Abstract 

en ' The Coulomb gauge model of QCD is studied with the introduction of 

<N : 

^^ . a confining potential into the scalar part of the vector potential. Using 

Q\ • a Green function formalism, we derive the self-energy for this model, 

i-G I which has both scalar and vector parts, S'^(p) and S^(p). A rotation 

Oh! 

JL . of these variables leads to the so-called gap and energy equations. We 

Mh. 

^ ■ then analyse the divergence structure of these equations. As this de- 

^ I pends explicitly on the form of potential, we give as examples both 

^ . the linear plus Coulomb and quadratically confining potentials. The 

C^ ' nature of the confining single particle Green function is investigated, 

and shown to be divergent due to the infrared singularities caused by 

the confining potential. Solutions to the gap equation for the simpler 

case of quadratic confinement are found both semi-analytically and 

numerically. At finite temperatures, the coupled set of equations are 

solved numerically in two decoupling approximations. Although chiral 

symmetry is found only to be exactly restored as T — > cxd, the chiral 

condensate displays a steep drop over a somewhat small temperature 

range. 



I. INTRODUCTION 

From a phenomeno logical and experimental viewpoint, two aspects of quantum 
chromodynamics (QCD) are believed to be fundamental, (a) that cliiral symmetry 
is spontaneously broken in the ground state and (b) that partons are confined. Both 
chiral symmetry and the phenomenon of confinement underly a common specula- 
tion, viz. that each is believed to give rise to a phase transition at a characteristic 
temperature or density T^{Tc) and PxiPc) respectively. As a side remark, one may 
note that these two phase transitions have a very different character - the chiral 
order parameter is zero in the ordered phase, and non-zero in the broken phase, 
in character with the notions of condensed matter physics, while the confinement 
transition has an inverted structure: the non-zero value of the order parameter, or 
Polyakov line in this case, usually characterizes the deconfined or ordered phase. 
The current belief, stemming from lattice gauge simulations, is that, for QCD in- 
cluding quarks, the critical temperatures coincide |1[] . This feature is not understood 
on a firm basis, and is merely an empirical observation. 

Phenomenological studies over the last decade have predominantly investigated 
the chiral symmetry aspect of QCD. For example, the Nambu-Jona-Lasinio (NJL) 
model P,^ has enjoyed tremendous success, due to its mathematical tractability and 
ease of use. Calculations of the low-energy mesonic and baryonic sectors demonstrate 
the well-known fact that chiral symmetry is broken at T = in nature. From a the- 
oretical point of view, one can study the model at finite temperatures and densities, 
but it is clear that the thermodynamical potential and associated bulk quantities 
will be dominated by inadmissable quark degrees of freedom in an intermediate tem- 
perature region where they should not exist at all ^. One conceivable method of 
removing this problem is to start with a confining model in the first place. This in 
itself presents many new difficulties that are not inherent in the non-confining NJL 
model. The fundamental theoretical questions that must be answered are (i) what is 
the analytic structure of a Green function that describes a confined particle ||^ , and 
(ii) what is the role and connection of chiral symmetry, if any, with confinement and 
the deconfining transition. In addition, such a model must be required to reproduce 
the excellent agreement of non-confining chiral models with experiment, which in 
itself, is no mean task due to the added complexity. Yet a further complication 



is that the nature of an effective deconfinenient potential is not precisely known. 
Bag model calculations suggest that confinement should be introduced in a scalar 
potential 0, while Coulomb gauge QCD naturally favors vector confinement. More 
detailed study is required to clarify this issue fully. Finally, even the understanding 
of these features still leaves the question open as to how the model should c?econfine 
dynamically. 

Several authors have attempted to address the problem that is associated with 
(i) above, i.e. of the analytic structure of a Green function for a confined quark [^0. 
Because of the various problems that are connected with these attempts (e.g. non- 
causality in the case of 0], and so on), this question is still not completely settled. 
In this paper, we do not make any assumptions about the form of the single quark 
Green function. Rather, we choose to study a model Lagrangian that contains a con- 
fining potential, and directly compute the Green function, assuming that it satisfies 
a Dyson-like equation. The feature of confinement should reflect itself automatically 
in the Green function and our task is to make this evident. Choosing a form for a 
confining interaction is not unique however. We have chosen to study the Coulomb 
gauge model that employs vector confinement as a starting point, simply because 
this model has been studied by several authors previously [P-P^. As such, this work 
is not entirely new: however, it unites previously unrelated work, and presents this 
model in a conceptual fashion that is different from most interpretations. Differences 
in the existing literature are clarified. The temperature dependence in this model is 
also examined. Here there are two different approaches that are possible, either that 
of Green functions in the real or imaginary time formalism, see for example [O], or 



via a variational principle that is applied to the free energy ||T^ . We use the former 
approach - one can explicitly see that it leads to equations that are entirely equiva- 
lent to those obtained in the latter approach. Here there has been some conflict in 
the literature. The gap equation that we have derived in this fashion corresponds to 



that derived by Ref. ||I3 on varying the free energy, but is in conflict with that of Ref. 



TT[| who disregard the temperature dependence of the quasiparticle energy, as well 



as that of Ref. [|1^ who simply write down a different gap equation. Further conflict 



arises also in the interpretation of the quasiparticle energy. In this paper, we utilize 
the knowledge gained from the invariance structure at zero temperature, T = 0, and 
implement a similar procedure at finite temperatures. Central to our approach is 



the appearance of physical quasiparticle energies that are temperature dependent, 
as is the case in the BCS theory of superconductivity for low temperatures. It is 
essentially this temperature dependence which is capable of driving a phase tran- 
sition, and which is studied here. This is in apposition to models that introduce 
an effective temperature dependence in the coupling strength, see, for example Ref. 



T3[| , which is difficult to deal with consistently from a thermodynamical point of 
view. In the case studied here, we keep the coupling constant fixed and examine 
the consequences of the temperature dependent quasiparticle energies. Numerical 
solutions are given for the coupled set of equations that form the gap equation in 
this case. It is interesting to note that although in part only approximate solutions 
can be given, the model dependence of {ipip) closely resembles that of the lattice 



simulations pO|]- Chiral symmetry, is however, not restored exactly. It would thus 
be of great interest in the future to have lattice simulations that could quantify the 
range of the drop in temperature of the chiral order parameter. 

Finally, one should note that the Coulomb gauge model, while appearing in- 
tuitively physical, is by definition no longer manifestly gauge covariant, as is the 
non-confining NJL model. However, since any potential introduced into the NJL 
model would also be expected to be instantaneous, this covariance is unavoidably 
lost. Thus a study of the Coulomb gauge model should provide many features that 
should be generic to four-fermion interactions and should be similarly handled in, 
for example, an extended NJL model that includes confinement. 

We point out that currently there has been a surge of interest in the problem 
of confinement. For different approaches, we refer the interested reader also to the 
works of Refs. [J18|,|19|,^T]-|23| . This list of references is, however, not exhaustive. 

This paper is organized as follows. In Section II, the general formalism is devel- 
oped for T = for an arbitrary potential. The form of the single particle propagator 
is given and the distinction between observables and non-observables is made accord- 
ing to the conditions of ultraviolet renormalizability as well as infrared finiteness. 
These divergences are discussed for two particular types of potentials, the Coulomb 
plus linear rising potential and the quadratically confining potential. In Section III, 
the solutions of the gap equation for a quadratically confining potential are discussed 
in this framework. In Section IV, the extension to finite temperatures is given, and 
numerical solutions for the condensate and dynamical mass is given. We summarize 



and conclude in Section V. 

II. FORMALISM FOR T = 
A. Self-consistent self-energy 

Our starting point is the Coulomb gauge QCD Hamiltonian, given by 

d^x'tp\x){-id ■ V)'^{x) (2.1) 



where Xa are the standard Gell-Mann matrices, and V{x) is an arbitrary potential 
at this point, and which will be considered to be confining. No current quark masses 
are considered. 

In general, the self-consistent self-energy S(p) associated with the interaction is a 
function of momentum. Under the assumption that the quark propagator is diagonal 
and independent of color, S{k)c,c' = S{k)6c,c', one has the mathematical consequence 
that the Hartree term vanishes. Thus the leading non-vanishing contribution to the 
self-energy, calculated to lowest order in the interaction strength follows from the 
Fock diagram that is depicted in Fig.l. One finds 



T.jT^.ViP-'kh^^S(ky,„^. (2.2) 



where S{k) is the self-consistently evaluated one particle fermion Green function, of 
form 

S{p)-'=^-E{p). (2.3) 

Note that this structure is particular to this form of the interaction, and, as such, 
may not be taken as being generic to confinement per se. Lorentz symmetry and 
invariance under a parity transformation enable one to decompose S(p) into a scalar 
and vector term, 

S(p-) = S^(b1)+7-pS^(b1), (2.4) 



where p = f>/\'p\ is the unit vector. Inserting the ansatz Eq. (^I^) into Eq. (|2.3|) leads 
to the exphcit form 

b[PO,P)- p2_[(^+sV^(p))2 + (s5(p))2] ^^■^) 



for the propagator, which, when inserted into Eq.( |2.2| ), leads to the coupled set of 
non-linear integral equations, 



^s,^ 2 f d^k ^^^^ ^^ S^(A:) 



^ ^^^~ 3J (27r)3^^^ ^^ p + S^(A;))2 + (S^(A;))2]i/2 ^^'^^ 

^ ^^^~ 3J (27r)3 ^ ^P "> [{k + ^y{k)Y + (S^(/t))2]i/2 ^^- ^^ 

for the unknown scalar and vector self-energies. The form of Eq.( |2.5| ) and also ( ^.71) 
suggests a reparametrization of S'^(p) and TX {p) as 

T.^{p)=E{p)sin<j){p) (2.8a) 

and p + TX (yp) = E{p) cos0(p) (2.8b) 

introducing the variable E{p) and the angular variable 0(p). It can be seen that the 
variable E{p) plays the role of an energy, by writing the propagator with respect to 
the new variables: 

7% - 7 ■ pE{p) cos{(f){p)) + E{p) sin(0(p)) 
Po - E{pY 



while the Eqs. (p.8|) are replaced by 

£;(p) sin(</>(p)) = -j j^V{p- k) sin{ct^{k)) (2.10a) 

E(p) cos(0(p)) -P=^j 7^^(P - ^) cos(0(A:))p ■ k. (2.10b) 

A simple manipulation of Eqs. (|2.10|) can be made to decouple the unknowns E{p) 
and (pip), leading to the result 

2 r d^k ^ 

psin(0(p)) = -j — -— l^(p- A;)[sin(0(A;))cos(0(p)) - p ■ fccos(0(A;)) sin(0(p))], 

(2.11) 

that determines 0(p) solely. One may note that the chiral condensate can be ex- 
pressed purely in terms of 0(p) as 

6 



(^^) = -A/ dkk''sm{^{k)), (2.12) 



vr^ Jo 



from which it is clear that 0(p) is the order parameter characterizing a non-trivial 
ground state. The determining equation for this function, Eq.( |2.11| ) is thus the 
gap equation for an arbitrary potential in this model. Assuming that the symmetry 
breaking is maximal for the particle at rest, one may impose the boundary conditions 



while for large p, 



lim0(p) = ^, (2.13) 



lim d)(p) -^ 0. (2.14) 



This completely specifies the problem. 

B. Confinement, physical and unphysical variables 

While the formal manipulations of the previous section are simple, some inter- 
pretation of the derived equations is required for a potential that is confining. In 
particular, the simple form of the propagator in Eq.( p.9|) might lead one to believe 
that this function has singularities in the complex plane. In fact, this is not so. 
Closer inspection of the defining relation for E{p) in Eq.( p.lO| ) indicate that it is 
related to the integral of the function V{p — k), which, for a confining potential, 
is infrared singular. Thus E{p) may not be interpreted as the physical energy of 
a quark. Only variables that are found to be infrared finite (and, if necessary, 
ultraviolet renormalizable) can be regarded in a physical sense ||10|| . 

We thus assert that: (i) E{p) is, in general, divergent, and thus necessitates 
a redefinition of the physical energy in order to extract a physical quasiparticle 
mass, and (ii) (j){p) is a physical quantity for both linear and quadratic confining 
potentials, as the gap equation that it satisfies is infrared finite. In the following two 
subsections, we examine the possible divergences that may be associated with the 
infrared, k ^ and ultraviolet A; — > oo behavior of the gap and energy equations 
respectively. 



1. Ultraviolet divergences 

a. Gap Equation For p — > oo, sin(0(p)) — > and cos(0(p)) -^ 1. The first term 
of the gap equation, Eq.( p.ll| ) is thus ultraviolet convergent. A divergence could 
however arise from the second term, depending on the form of the potential. For a 
Coulomb-like potential of the form V{r) ~ — 1/r or V{q) ~ l/q^, a divergence going 
as 

d'k^,-jdke^^ (2.15) 

occurs. In general, the gap equation is ultraviolet divergent for potentials that fall 
off more weakly than the Coulomb potential. In order to quantify this statement, 
consider a potential of the general form V = V{\p — A;p). The leading k dependence 
is obtained from the Taylor series of this function about p = 0, as 

V{\p- k\') = V{k') - 2p. k^^01 + ... (2.16) 

and inserting this into the divergent term of the gap equation 

^J j^,V{\p-k\')p.ksin{<p{p)), (2.17) 

one sees, after performing the angular integration that the leading possibly divergent 
contribution goes as 

For a confining potential that has the limiting behavior 

y(g2) ^~°^ (g^)-", a>0, (2.19) 



Eq.( p^.l8[ ) can be evaluated to give the result k ^""'"^/(— 2a + 2), which thus exhibits 



a divergence for a < 1. Denoting 

^ ^ A f d^k dv{e] 



k, (2.20) 



9 7 (27r)3 d{k^] 
one finds that the divergent part of the gap equation, i.e. Eq.( |2.17D , is of the form 



divergence = {Z — l)psin(0(p)). (2.21) 
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This, however, has exactly the form that appears on the left hand side of the gap 
equation, and thus any ultraviolet divergences can be simply compensated by intro- 
ducing an appropriate counter term into the Hamiltonian, in this case, 

AH = (Z-l) I d^xi)\x){-ia-V)i){x). (2.22) 

With this renormalization prescription, ultraviolet divergences can be unambigu- 
ously dealt with. 

h. Energy equation Since 0(A;) -^ when k -^ oo, it follows from Eq.( p.lOD 
that Eij)) is ultraviolet convergent. 

2. Infrared divergences 

a. Gap equation By performing an expansion of the bracketed term in the 
integrand of Eq.( |2.11| ) about p = k, i.e. 



sin(0(fc)) cos{(f){p)) — p ■ kcos{(f>{k)) sin(0(p)) 

= (P - k) ■ p<P\p) + 0((p - kf), (2.23) 

and shifting the variable q = p — k in the integrand of the gap equation one may 
isolate the possible divergence as being 



divergence ^^J Al_V{q)0{q^). (2.24) 

Thus, for a linear confining potential, V{r) ~ r or V{q) ~ 1/q'^, this term takes the 
form 

/ d\-^0{q') ~ j dqq'\0iq') ~ J dq, (2.25) 

which is infrared finite. Thus the gap equation is non-singular in the infrared in 
this case. This can also be verified to be true for a potential that is quadratically 
confining, i.e. V{r) ~ r^. 

b. Energy equation Eq. (|2.10|) can be used to determine the quasiparticle energy. 



2 f d^ ^ sin(0(fc)) 2 f d'^k ,,,,^, sm{(f){\p-k\)) 

^^^) = 37 (2^^^"^- ^");R0M) = 37 (2^^^^) sin(0(p)) " ^2-26) 



The divergence structure in the infrared follows on expanding the sine function about 
k = p. The leading term is 



2 r d^k 
divergent part of i?(p) ~ - / - — —Vik) 

3 J (27r)'^ 



d?k 
divergent part of E(0), (2.27) 



the latter statement of which can be seen by inspection. It is thus evident that the 
energy difference 

E{p) - E{0) (2.28) 

is infrared finite. 

One may note that the subtraction of the infinite zero point energy corresponds 
to shifting the potential by an infinite constant 

2 r d^k 
nrl^V{r)~-J j^^Vik) (2.29) 

or 

V{q) - V{q) - i'^^f^^'\^lj^y(k)- (2-30) 

Thus the content of this section on infrared divergences can also be viewed in the 
context that physical observables are invariant under shifts of the potential by a 
constant. This is immediately obvious for the gap equation, Eq.( |2.11| ), but not for 
E{p) alone as in Eq.( ^.2(j| ). Under the general shift 

V{f) -^ V{r) + [/, (2.31) 

it follows that E{p) -^ E{p) + 2f//3. Consequently, the determination of -E'(O) does 
not provide a unique determination of the quasiparticle mass. Following Ref. ||10|| , 
one may make an identification of the quasiparticle mass by rewriting the single 
particle Green function of Eq.( |2.5|) as 



SU, ^ = ?±M . ^zM (2 32) 

^^°'^^ Po + iEip)-EiO)) + EiO))^Po-iEip)~EiO))-EiO)y ^^ ^ 



where the residues 



R±{p) = \{loTl ■ pcos(0(p)) ± sin(0(p))) (2.33) 
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are unaffected by a shift in the potential. By comparing the residues at low momenta 
with those of a free relativistic particle of mass m*, viz. R\ = ^(7° =p 7 ■ p/m* ± 1), 
one arrives at an unambiguous identification 

J_ ^ ^^ cos{<P{p)) ^2.34) 

m* p^o p 

for the quasiparticle mass. 

It is however physically more intuitive to introduce the shifted energy 

E{p) = E{p) - E(0) + m* (2.35) 

which has the desirable property -E'(O) = m*. In order to ensure this property, the 
shift 

-U = -E{0)+m* (2.36) 

must be invoked in the potential Eq.( p.31| ). 

In concluding this section, we point out that both m* and the condensate (ipip), 
calculated via Eqs.( p.34 ) and ( p. 12 ) respectively, fullfil the requirements of being 



physical observables as they are solely functions of the gap variable 0. 

III. SOLUTIONS FOR QUADRATIC CONFINEMENT AT T = 

In this section, we illustrate the solution of the gap equation for a quadratically 
confining potential, since such a form for the potential has the advantage that the 
integral gap and energy equations can be reduced to differential forms, even although 
they are highly non-linear, and the problem is a boundary value, not an initial 
value one. Nevertheless, many analytic techniques are at ones disposal, and in 
this case, precise solutions can be found. Note that this problem has been solved 



approximately elsewhere [16 , but not within the physical framework presented here. 



Here we give a more general solution than that in Ref. |]I6|. Since the method will 
also be applied to the problem at finite temperatures in Section IV, we do this in 
some detail. 

Our starting point is the attractive potential 

V{f) = -Xr^ + U, (3.1) 

11 



with A > controlling the effective strength of the potential, and U determined via 
the condition Eq.( |2.36|) . This potential is binding in the quark- ant iquark channel. 
The Fourier transform of the potential Eq.( p.l|) is given as 

V{q) = \i27rf \/l5^^\q) + U {2-Kf 5^^\q) . (3.2) 

Inserting this into the integral gap equation, Eq.( p.ll[) and performing an integration 
by parts, leads one to the differential form 

psin(0(p)) = \\{(t>"{p) + -0'(p) + 4 sm(20(p))), (3.3) 

3 p p'^ 



which is to be solved subject to the boundary conditions Eqs. (|2.13|) and ( p.l4| ). A 



semi-analytic solution to this problem is presented in the following subsection. This 
is later compared with an exact numerical solution. 

Semi-analytic solution 

The differential gap equation, written as 

ap^ sin(0(p)) = p'^(j)"{p) + 2p(f)'{p) + sin(20(p)) (3.4) 

with a~^ = 4A/3 > can be solved by the technique of asymptotic matching p4| . 



In this method, one identifies overlapping regions in the independent variable, for 
which exact analytic solutions can usually be obtained. On enforcing the boundary 
conditions, these solutions are matched onto one another asymptotically in their 
common region of validity. Four such regions can be identified, of which the second 
and fourth are a limiting cases of the third and are thus not independent. These 
regions are 

(1) p^O =^ p^(t)'\p) + 2p<p\p) = - sin(20(p)) (3.5a) 

(2) p and (/)(p) small =^ P^4'"{p) + '^P4'\p) = ~'^4'{p) (3.5b) 

(3) (j){p) small for arbitrary p =^ p'^4>"{p) + '^P4''{.P) = ('^P^ ~ 2)(/)(p) (3.5c) 

(4) (j){p) small and p — » oo => P^<P"{p) + '^p4'\p) = ap^4>{p)- (3.5d) 

Correspondingly, the energy equation, in differential form, is given as 

E(p) = --(0'(p)2 + lcos2(0(p)) +pcos(0(p)) + ^U. (3.6) 

a p^ 3 

Thus, once one has knowledge of 0(p), the energy may be immediately determined. 

We now deal with each of the above cases in turn. 
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case 1 : The determining equation Eq. (|3.5| a) in this case cannot be solved 
exactly due to non-linearities. However, since the boundary condition Eq.( |2.13D 



must be imposed, it is useful to write 0(p) = 7r/2 + e{p) with e(0) = 0. One may 
insert this ansatz into Eq.( p.5| a) and solve the resulting equation exactly for small 
e. The non-divergent solution is found to be linear in p so that one may write, in 
leading order, 

(t>{p)^\ + cp + r,{j>), p^O, r/(0) = (3.7) 

where c is as yet an undetermined constant. Clearly this procedure can be repeated 
by constructing the differential equation for 77 (p) to obtain further corrections. This 
leads to the higher order correction ri{p) = bp^ + 0{p'^), with b = {a — 4c^/3)/10. 
Note that the constant c in Eq.( |3.7|) plays an extremely important role. According 
toEq.(^) andEq.Q, 

m*-^ = -(/)'(0) = -c, (3.8) 

i.e. the initial slope of the gap angle is related to the inverse quasiparticle mass. 
This determines c to be negative, c = — |c|. 

We now turn to the determination of E{p). Using the approximation Eq.( |3.7| ) in 
Eq.( p.6| ), one finds 

E(p) = -— + [—- 2c]p2 + |t/, p^O. (3.9) 

d CI o 

The criterion that E{0) = m* now fixes U to be 

C/ = ^^ + 51. (3.10) 

It is interesting to note that, in light of these comments, the dispersion relation for 
E{p) for small momenta, from Eq.(3.8), 

E{p)=m* + Bp'^, p^O, (3.11) 

with B = 2c'*/a — 2c > 0, is non-relativistic. 

case 2 : In the region specified as (2), the determining equation Eq.( |3.5| b) can 
be solved exactly. The solution is 

^,(p) = J^cos{^\n{p) + 6) (3.12) 

where 7 and 6 are two constants of integration. 

13 



^s{p)=Bp-'/'K^^^,{-V^p'/'), (3.13) 



case 3 : The determining equation over this region, i.e. Eq.( p.5| c) is a Bessel 
equation. Selecting the solution that falls off with momentum leads to the solution 

,2 
3 

where Kn{x) is the nth order modified Bessel function of the second kind, and B is 
an integration constant. 

case 4 '■ Eq (|3.5| d) is a limiting case of Eq. (|3.5| c) and has as decaying solution, 
the modified Bessel function 

Mp) = Cp-'^'K^/,i^V^p'/'). (3.14) 

For large x, Kn{x) ~ xf^^^^V- + 0{l/x)], so that 

04(p) ~ Ap-''/^ exp(-^v^p3/2) ^3^^5) 

in the limit p ^ oo. 

We now turn to the energy determination. From Eq.( |3.6| ), it follows that 

E{p)''^"^p. (3.16) 

Thus one recovers a relativistic dispersion relation at high momenta. 

The next issue is to match all solutions that have been obtained piecewise to 
construct one continuous solution over the entire momentum range. 

Asymptotic matching over the intervals. One may observe that the Eqs.( ^.5| a) 
and (b) are scale invariant, i.e. invariant under the transformation p — > \c\p. This 
enables one to determine the constants 7 and 5 for the solution (|3.12| ) without 
actually knowing c = </>'(0). One determines the solution to Eq.( |3.5| a) with initial 
values 

</>i(0) = ^ and 0/(0) = -1, (3.17) 

and one considers this solution for p -^ 00. This is matched onto the solution 
02(p) from Eq.( |3.12|) and the associated 4>'2{p) in the regime p ^ 00 of the phii{p) 
equation. This is possible since the solution of 02 (p) can be matched to that of 0i(j5) 
down to small momenta. On the other hand, both Eqs.( |3.5| a) and ( |3.5| b) have the 
same behavior in the limit p -^ 0. One is thus able to extract the constants for the 
solution in region (2). 

14 



Numerically, one finds 

7 = 0.700460a-^/'' 

5 = 0.238503. (3.18) 

Now the solution in region (2), 4'2{p), niust be matched onto the solution in region 
(3), 03 (p). To do this, one requires 03 (p), for p ^ 0. Since 

iri,"~°ir(zz/)(ia;)-^^ + c.c., (3.19) 

one finds 

Mp) '- Bp-'/Yi^)\l[eM^^ H^P'^') - ^P) + cc] (3.20) 

= 5p-^/^|r(^)|cos(^ln(p) + ^ln(^)-p), (3.21) 

where p = arg(r(^)) ~ —1.877978. The full, non-scaled solution in region (2), 

lip) = -^ cos(^ ln(|cb) + 6) (3.22) 

\c\p ^ 



can be directly compared with Eq.( |3.2lD , to give 



and 



5 + ^ In |c| = ^ In(^) - p + nvr (3.23) 



^ =i?|r(z^)|. (3.24) 



One may now extract the slope \c\ from Eq. (|3.23 ). Explicitly, one has 



a^/3 2 

|c| = -^exp^(r27r-p-(5). (3.25) 

where n is the number of nodes occurring in the solution. In Table I, we list the node 
and the value for c that were extracted via this asymptotic matching procedure via 
expression ( p.25 ) in comparison with a trial exact procedure in which one attempts 



to solve the full equation starting as near to the origin as is numerically permissible. 

In Fig. 2, we show the nodeless solution for the gap parameter. The analytic 

solution in regions (1), (2) and (3) are indicated on the graph by various symbols. One 
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notes that the approximation is in excellent agreement with the numerical result. 
This procedure has proven extremely useful both as a check but also in guiding the 
numerical solutions, since the gap equation is highly non-linear. Figure 3 indicates 
the first three solutions for (f){p). One notes that the slope at the origin increases 
sharply with the number of nodes, indicating that a succession of values of m* exist 
that decrease with the number of nodes. 

Energy Density. In order to determine which of the solutions, i.e. that with 
no node present or alternatively with multiple nodes present should give rise to the 
ground state of the system, one needs to evaluate the change in the energy density 
that occurs in going from the normal vacuum to the condensed one. The physical 
solution to the gap equation is then determined as the solution that lowers the 
energy the most. The vacuum energy density in the condensed phase, W condensed = 
{cond\H\cond) or in the normal phase, Wnormai = (0|-ff|0), can be obtained directly 
from the appropriate single particle Green's function, via the formula [^] 

W^cn = -^ / ^^^[7% + 7 ■ Pl5c,n(Po,p)e*''^° (3.26) 

where c, n refer to the condensed or normal phases respectively. This equation is 
general and exact for a four-point or two-body interaction, such as that under study 
here. Explicitly inserting the Green function, Eq.( p.32| ), for the condensed phase 
and the free Green function for the normal phase, one finds the expression 



Q /-oo 11 

^W = Wc-Wr, = —^ dpp'[2p{l - cos(0(p))) - — sin2 0(p) + — 
ZTT^ JO ap^ 2a 



(3.27) 



for the change in energy density in moving from the normal to condensed vacua. A 
numerical analysis indicates that ISW is negative, and that the nodeless solution, 
i.e. the solution with the smallest slope or largest mass, always leads to the lowest 
energy state. 

IV. FINITE TEMPERATURE 

To date, there have been several studies of the finite temperature gap equation. 
These differ markedly in (a) the gap equation actually proposed or derived and the 
subsequent interpretation of the variables, and (b) the form of the potential chosen 
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for studying this model. In the first part (a), there can of course be only one correct 
equation. In analogy to the T = case, it can be formally derived by one of two 
equivalent methods, the Green function approach, similar to that followed at T = 
in Section II of this paper, or via a variational method that proceeds by constructing 
the Gibbs free energy and minimizing it with respect to a variational parameter, see 



Ref. [0. In this section, we will take the former approach, since the generalization 



from Section II is evidently simple. We contrast this result with the conflicting 



postulated gap equations of other authors [0 and interpret the quasiparticle energy 



that occurs along the lines of Refs. |T^ and |T^|. 

With regard to point (b) above, let us mention that the only attempt to date 
of extracting a numerical solution to a finite temperature gap equation (not that 
derived here) employs an a priori temperature dependent coupling strength for a 
Coulomb-like potential |]13[. This, in itself, presents difficulties in the construction of 
a consistent thermodynamic description. In the case handled, the second order phase 
transition that was observed was due to the coupling, that strongly modifies the 
infrared behavior of the potential, rather than the inherent temperature dependence 
of the quasiparticle energy, which was omitted in this description. 

An additional difficulty manifests itself in the finite temperature formalism with 
regard to the zero temperature one. Since statistical averages over a grand canon- 
ical ensemble are made, all states are included in the averaging process. However, 
physical states should be colorless. Imposing this restriction leads to rather complex 
technicalities and a complex set of equations, from which it is not evident what the 
physics actually is [|r^. Taking the point of view that the restriction to color singlet 
space should cause changes of at most 10 — 20%, but should not induce radical differ- 



ences otherwise [^, we attempt to solve the gap equations that are derived without 
this restriction and which have a physical constraint imposed upon the quasiparticle 
energy. We do this using the confining potential that was introduced in Section III, 
following the analysis of this section rather closely. 

We proceed by briefiy deriving the finite temperature gap equation. 
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A. Gap equation 

The Green function approach that was detailed in Section II can be simply 
generalized to finite temperatures in order to derive the associated gap equation. 
Here, using the Matsubara or imaginary time propagator, 

where Un = (2n + l)7r//3, n = 0, ±1, ±2, ±3, . . ., and the self-energy S^(A;) is again 
decomposable into a vector and scalar piece, S^(A;) = S^(A;) + 7 ■ A;S^(/i;), one may 
write 



18 „ jSl \a \a 

^M = —qT.T.J 7^3^(p~ kh\S,{^u^, fc)7oy , (4.2) 



on interpreting the Fock diagram of Fig. 2. Now the steps that were performed in 
Section II leading to Eq.( p.l(]| ) can be similarly executed here, once the Matsubara 
sum has been performed. Eqs. (|2.1CI| a) and ( p.lO| b) are replaced by their temperature 
generalizations, 

Ep{p) sin(0^(p)) = -j j^,V{p- k) sin(0^(A;)) tanh(|E^(A;)) (4.3a) 

E^ip) cos{Mp)) -P = ^J J^s^^P- ^)^ ■ ^ cos(0^(A;)) tanh(|E/,(A;)). (4.3b) 

Note that a central issue in this derivation is that the self-energies S^(p) and S^(p) 
making up "^fsip) are now necessarily a function of temperature. This temperature 
dependence is now relayed to the reparametrization functions 0/3 (p) and Ep{p) that 
are defined as in Eqs.( p.8|) , but which contain the explicit temperature labels. It has 
been noted that the simple procedure used previously to decouple these equations 
no longer functions |T^|T^. The finite temperature gap equation, contructed by 
multiplying Eq. (^l3| a) by cos(0/3(p)), Eq. (^l3| b) by sin(0/3(p)) and subtracting, is 

2 f d^k 
psm{(j)f3{p)) = -J -——V{p - k)[sm{(j)f^{k)) cos{(l)fs{p)) - p ■ kcos{(f)fs{k))sm{(j)p{p))] 

X tanh{(^Ep{k)). (4.4) 



This gap equation coincides with that of Refs. [|12| and [T^ . It differs however from 



that of Ref. |TT| in that these authors regard E^^ = E as a constant, independent of 
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the temperature. In contrast, in Ref. uM, the authors use a gap equation similar 



to Eg. ([4 .41) , but with the tanh factor entirely absent. These authors claim that the 
hyperbolic tangent function and energy dependence is irrelevant, arguing that E{k) 
is infinite anyway. This argument is, however, too naive, and is not confirmed by 
our study, in which we first make use of physical arguments to shift the energy along 
the lines that were discussed in Section II: Note that, taken at face value, Eq.( [4.4|) 
taken together with Eq. ([4.3| a), violates the criterion that physical observables, such 
as (p/sik) are invariant under a shift in the potential by a constant, as was easily seen 
to be the case in Section II at T = 0. This issue has been dealt with in detail in Refs. 
[12[| and [|17|. Both of these sets of authors show that this conflict can be resolved 



by supressing color fluctuations and restricting the grand canonical ensemble trace 
to cover color singlet states only. The color projected gap equation [|1^ is however 



extremely difficult to solve, in particular, because of the temperature dependence of 
the energy. Note that the temperature dependence of the energy variable is crucial in 
the BCS theory of superconductivity for obtaining the correct temperature behavior 
of the gap parameter for a superconductor [^|. It is thus the purpose of this study 



to investigate the role of the temperature dependence of the energy function and its 
consequences. 

Guided by the situation at T = 0, one may make a simple physical ansatz, that 
of replacing Epip) in the tanh function of Eq.( [4.4|) by 

Ep{p) = Ep{p)-Ep{Q)+m*^. (4.5) 

This is motivated, as in Section II following Ref. |T^, by comparing the Green 
function in this case with the low momentum version of the free temperature Green 



function 121. The real time forms are 



^(p , c.) = [^%M_ + R-^^^P) 



uj - Ep{p) +17] uj + Ep{p) -if] 
+ 2mf[l3Ep{p)] [70,; - 7 . pEp{p) cos(0^(p)) + Ep{p) sin(</)^(p))]5(^2 _ ^2^^)^^ 

(4.6) 
for the full Green function, with 

R±Ap) = -[7°T7-Pcos(0^(p)) ±sin(0/3(p))], (4.7) 
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where f{x) = l/{exp{j3x) + 1) and t] -^ 0+, while the free Green function has the 
form 

R%ip) , R%ip) 
uj — m*n + 11] u + ml — iTj 
+ 27rz/(/?m;)[7ocu - 7 ■ P + m*f,]5{u^ - m}) (4.8) 

with R%{j)) = 1(70 ^ 7 . p/m*p ± 1). 

Thus one recovers the identifying criterion 

hm^^^^M^ = ^ (4.9) 



that generahzes Eg. ( p. 341) to finite temperatures, as well as the essential condition 



Ep{p) ^ Ep{p) ^ m% p^O (4.10) 

that has prompted the replacement Eq. (^4.5| ). Given Eq.( [4.5| ), one recovers the ex- 
pected physical quasiparticle properties. 

The associated finite temperature condensate density is determined to be 



3 |-oo g 

-^ / (iA;/c^sin(0/3(A;)) tanh(— . 
vr^ Jo 2 



{iPiIj)p = -j dkk^sm{(t)p{k))tiinh{^Ep{k)), (4.11) 



which will be evaluated in the following section for the potential for quadratic con- 
finement. 

B. Quadratic Confinement 

The purpose of this section is to examine the temperature dependence of the gap 
function and the self-energies for the quadratically confining potential of Eq.( p.l| ) 
and its Fourier transform, Eq.( |3.2| ). In this section, the strategy is to implement 
Eq.( [4.5|) in an approximate fashion that highlights the temperature dependence of 
the quasiparticle energy. We proceed once again by bringing the gap equation into 
differential form. One finds that Eq.(p.4|) has the generalized form 



ap sm{(f)p{p)) coth{-Ep{p)) = p (p" p{j)) + 2p0^(p) + sin(20^(p)) 

+ ^P' ■ imp] vM "^- C-i^) 

smh{(3Ep{p)) 
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while the energy equation, Eq. (|3.6|) is generahzed to read 



a p^ 2 



Efiip) = --ih'ipf + -cos2(0^(p)))tanh(^E^(p)) +pcos(0^(p)) 



+ ^[pt^nh{^Ep{p))]" + ^f/tanh(^E^(p)). (4.13) 



Unhke the system of equations ( p.4| ) and ( p.6| ) at zero temperature, one sees that 
Eqs.( |4.1^ ) and ( [4.13D are coupled non-linear differential equations in the two vari- 



ables (j)j3 and Ej3. While it is possible to make a simple analysis for high temperatures 
and low/high momenta, it is a difficult problem to solve these two equations exactly. 
We thus do not attempt to do so here. After discussing the limits just mentioned, we 
shall investigate the role of an assumed temperature dependent dispersion relation 
for the energy, that allows us to decouple these equations and solve them. 

c. High temperature limit. If, as T ^ cxd, EjbIp) remains finite, then from 



Eq.(M2D 



psin(0^_^o(p)) = for all p (4-14) 

so that 

0/3^o(p) = 0. (4.15) 

Consequently, from Eq.( |4.13| ), it follows that Ep_^o = p for all p, and one recovers 
the chiral limit. 

d. Small momenta. Using the boundary conditions, Ei3{0) = mp and 0/3(0) = 
71/2, one can expand 0/3 (p) and Ef^{p) in the small momentum region, 

0/3 (P) = ^ + ci P + (ip^ + ■■■ (4- 16a) 

E(i{p) = mp + a^p + a^p'^ + ... (4.16b) 

where, from Eq.( |4.1J| ), 

4 = -^ = C/3. (4.17) 



One may insert these expansions into Eq.( [4.13D . A necessary condition that Ep{p) 
does not diverge as p ^ is that a^ = 0. Thus one regains the observed T — > 
dependence of Section III that 
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Ef^ip) = m*f^ + a^p"^, for p ^ 0, (4.18) 

this equation being the analog of Eq. p.ll| ). As in that case, one may continue the 
expansion for 0/3 (p), writing 

<P(s{p) = \ + c^P + Mp), with r]p{0) = 0. (4.19) 

Inserting this and Eq. ([4.16| b) into the gap equation enables one to find a differential 
equation for 77/3 (p) which has the solution rjjs^p) = b^p^, with 

^^l|.cotlW2|c^|,.-iggl^^ 1,1^1)3], (4.20) 

The coefficient 02 of the energy expansion may also be determined by inserting the 
expansions for E and into Eq.( |4.13| ). One finds 

4 = l[3|c'^|2tanh(/3/2|c^|) + a(^ - ^f/tanh(/5/2|c^|))]i cosh2(/3/2|c'3|). (4.21) 

Thus the behavior of (pjs and Ef^ for small momenta is similar to that obtained at 
T = 0, see Section III. 

e. Large momenta. From Eq.( [4.13| ), it follows that, at large momenta, Ep{p) ~ 
p, p -^ 00, and one no longer has a temperature dependence in this quantity. Thus 
the gap equation is precisely that of the zero temperature case. 

As has been seen at zero temperature, and now also at finite values of the temper- 
ature, the dispersion relation for the energy is non-relativistic in the low momentum 
regime, but relativistic for high momenta. This is commensurate with a general 
relativistic ansatz for the quasiparticle energy. In the following, we will discuss the 
solution to the gap equation for two possible ansdtze: (i) momentum independent, 
but temperature dependent, 

Mp) = m; (4.22) 

and (n) the relativistic dispersion relation. 



E^{p) = ^m*p^+p^. (4.23) 

Note that ansatz (i), being totally devoid of momentum dependence, serves to high- 
light the infiuence of the tem,perature dependence of the energy variable on deter- 
mining the features of the solution to the gap equation. In addition, the simplicity 
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of this ansatz allows us to solve the resulting equations easily. The second ansatz, 
while meeting physical expectations more readily, leads to equations that are some- 
what more difficult to solve. As shall become evident, the results obtained from the 
second ansatz do not differ qualitatively from that of ansatz (i). This appears to 
indicate that the momentum dependence is secondary in comparison to a correct 
treatment of the implicit temperature behavior of the quasiparticle energy. Our two 
choices are detailed in the following two subsections. 

C. Momentum independent ansatz 

In this subsection, we examine the solution of the temperature dependent gap 
equation under the ansatz Eq. (|4.22|) . This is interesting to study since it contains 



the essential ingredient, i.e. temperature dependence, that has previously been 
overlooked. Furthermore, the resulting transcendental equation is simple to solve. 
Let us be more precise. Inserting Eq.( ^.22D into Eq.( ^.12| ) leads to the result 



aV MMp)) = p'<P'/3'ip) + '^pMp) + sin(20^(p)), (4.24) 

with 

a^ = acoth(-m^). (4.25) 

Except for the cotangent factor, this equation is identical with the T = equation, 
Eq.( p.4|) . It may thus be treated semi-analytically in the same fashion. In particular, 
the relation Eq.( |3.25|) is valid, with c replaced by c^ and a replaced by a^, i.e. 



\c^\ = L^exp(-^(n7r -p-6)). (4.26) 



Now, inserting a^ from Eq. (^4.25|) into this equation leads to the simple result 

(coth(fm*))i/='|c|, (4.27) 



which, together with identification jc'^l = ^/Tn*g and \c\ = 1/m*, leads to the tran- 
scendental equation 



mj = m*(tanh(-m^))^/^ (4.28) 
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that determines ml. In addition to the trivial solution, ml = 0, this equation has 
a non-trivial solution that can be determined numerically. Following this, (ppip) is 
determined, and the condensate {ipip), from Eg. ( [4. Ill) , can be shown to have the 
form 

{i,i^)p = (^^)T=o(tanh(|m;))2 (4.29) 

in relation to its zero temperature value. 

In Fig. 4, we show the values of (pisip) as a function of p for different values 
of the temperature, scaled with respect to the dynamically generated constituent 
quark mass. One notes from this figure that the slope increases with increasing 
temperature, indicating the fall off of the constituent quark mass. The temperature 
dependence of the condensate is shown in Fig. 5 as a function of the dimensionless 
variable T/m*. One sees that it falls off strongly to zero as T — > 0: at T/m* = 2.0, 
the value of the condensate is less than 5% of its value at the origin. However, the 
transition of the condensate is smooth and there is no evident critical temperature 
at which a phase transition that restores chiral symmetry occurs. This behavior is 
also evident in that of the quasiparticle mass, which is shown in Fig. 5 as a function 
of T/m*. This function decreases rapidly over the first unit of T/m,*, but thereafter 
drifts gradually to zero. 

D. Relativistic dispersion relation 

In this subsection, we investigate the use of a relativistic energy- momentum 
ansatz, viz. that of Eq. ( [4.231) . For small values of p, this has the expansion 



Mp)'='m*, + ^ + 0{p'), (4.30) 

enabling one to identify the expansion coefficient of the energy as 

4 = J- = -\c% (4.31) 

^ 2m*p 2' '' ^ ^ 

and that of the gap parameter, from Eq.( |4.2U( ). One then has 
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The estimate for the gap parameter, (pjslp), for p small, thus enables one to perform 
a numerical integration starting at a value of p close to the origin. Fig. 6 shows the 
solution to the gap equation as a function of momentum at different values of the 
temperature. This bears close resemblance to the solution obtained in the previous 
case. The condensate, calculated numerically from Eq.( ^.ll| ), is shown in Fig. 7. As 
in the previous case, it goes to zero, but somewhat faster than in the previous case. 
Once again, no transition to a chirally symmetric phase is observed at a finite value 
of the temperature, although the fall off, for low temperatures, is steep. 

E. Comparative analysis 

It is useful to examine the results for the previous two calculations on one graph. 
In Fig. 8, we show the condensate as a function of the temperature, for two arbitrary 
choices of the mass at zero temperature, either m* = 200MeV, or m* = SOOMeV. 
It is interesting to note that both assumptions for the dispersion relation lead to 
the same qualitative picture. Both curves fall rapidly initially. Instead of a phase 
transition occurring however, the curves flatten off and fall more gently to zero. One 
may speculate that a better form or exact solution of the coupled equations may 
in fact yield a chiral phase transition at a finite value of the temperature. In any 
event, it is heuristically evident that the temperature fluctuations alone could drive 
a phase transition. In view of the difficulties presented here, we must conclude that 
it is necessary to study temperature behavior of even simpler confining models. 

The temperature behavior of the masses themselves is shown in Fig. 9. These 
are also decreasing with temperature, as occurs, for example, in the Nambu-Jona- 
Lasinio model, but no direct linear proportionality with the condensate is in ev- 
idence. The fact, however, that the quasiparticle masses are a function of tem- 
perature is an indication that this feature, which gives rise to dropping or rising 
composite particle masses, will continue to be an existing feature of a model that 
incorporates confinement and chiral symmetry and which one can hopefully solve 
exactly. It is to be believed that the solutions obtained here, which in both cases in- 
dicate a drop in the condensate energy, but no phase transition at finite temperature, 
are a consequence of the approximations that have been made. 
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V. SUMMARY AND CONCLUSIONS 

In this paper, we have analysed the Coulomb gauge model with respect to the 
inclusion of a divergent potential. In this model, confinement is incorporated into 
the scalar sector of the vector four-potential, although there is no compelling reason 
why this should be done. The consequences for the quark propagator are discussed 
and an analysis of the divergences in the infrared sector as well as the ultraviolet is 



presented. The arguments of Adler and Davis [|1^ that physical quantities should 
be infrared finite, coinciding with the condition that a shift in the potential by a 
constant should not affect physical quantities for color singlet states is used as the 
determining criterion for physical observables. An interpretation of the quasiparticle 
mass is made, since the quasiparticle energy is necessarily divergent. 

The finite temperature generalization of this model yields a gap equation which 
differs from that of several authors. The identification of the quasiparticle mass in 
the fashion that was done at T = leads to a natural physical prescription that 
suggests the introduction of a non-divergent temperature dependent quasiparticle 
energy Epij)) = E{p) — E{0) + m*^, so that Ep{j)) — > m^ as p — > 0. Even so, the 
complex set of coupled equations that one obtains for the gap parameter and the 
energy cannot be solved exactly. Since the issue of other authors has been to neglect 
the temperature dependence of the energy function and obtain alternatively that a 
phase transition can |T3] or cannot [^ occur, we propose to investigate two simple 
dispersion forms for the quasiparticle energy that have an inherent temperature 
dependence. We find that, while there is no transition to a chirally symmetric 
phase at a finite value of the temperature, there is nonetheless a sharp drop in the 
condensate density or order parameter at some point. Whether this is endemic to a 
confining potential, or whether it is simply a consequence of the shortcomings of the 
approximations made, is as yet unknown and could be studied using other confining 
potentials, and simplifying models. It is perhaps therefore of interest to develop 
models that are closer in structure to that of Nambu and Jona-Lasinio, in order to 
understand these properties. 
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TABLE I 



n 


1 Casmpt 


1 Cnurn \ 





2.09169 


2.03750 


1 


22.4837 


22.4246 


2 


241.680 


241.634 


3 


2597.84 


2597.83 


4 


27924.4 


27925.0 


5 


300161 


300170 



Table 1: Solutions for the slope c obtained via asymptotic matching of the ap- 
proximate solutions, and via a numerical procedure. The symbol n denotes the 
number of nodes that the solution has. 
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FIGURE CAPTIONS 

Fig. 1. Fock Feynman diagram for the self-energy. 

Fig. 2. The nodeless solution for the gap angular variable is shown as a function 
of the scaled momentum k = a^^^p. The semi-analytic solutions in regions (1), (2) 
and (3) are indicated by crosses, stars and circles respectively. 

Fig. 3. The first three solutions for (f){k) with zero, one and two nodes, are shown 
as a function of the scaled momentum k = a^^^p. 

Fig. 4. The quasiparticle spectrum Ej3{k) and the scalar and vector self-energies 
T.^{k) and S^(A;) + k are shown as a function of the scaled momentum k = a^^^p. 

Fig. 5. Solutions (ppip) as a function of p, in the approximation Ep = m*^ for 
various values of the scaled temperature T/m*. The upper curve corresponds to 
T/m* = 0. Successive curves have the values T/m* = 0.2, 0.5, 1.0, 5.0 and 10.0. 

Fig. 6. The temperature dependence of the condensate, normalized to its zero 
temperature value. The approximation Ep{p) = m*^ has been used. 

Fig. 7. The temperature dependence of the dynamically generated quark mass 
is given as a function of the scaled temperature T/m*. 

Fig. 8. As in Fig. 5, with the momentum dependent ansatz -E'|(p) = m,*^"^ + p^ . 

Fig. 9. As in Fig. 6, with the momentum dependent ansatz E'l{p) = m*^^ + p^ . 

Fig. 10. The condensate, shown for two arbitrary choices of the dynamically 
generated masses m* = 200MeV and m,* = 300MeV. The upper curves in each 
case correspond to the momentum independent ansatz Ep{p) = m,*^, while the lower 
correspond in each case to the dispersion ansatz E'i{p) = p^ + m*^^ 

Fig. 11. The current quark masses, shown for two arbitrary choices of the dy- 
namically generated masses m* = 200MeV and m,* = 300MeV. The upper curves in 
each case correspond to the momentum independent ansatz Ep[p) = m,*g, while the 
lower correspond to the dispersion ansatz E*Jp) = p^ + 1^*0^ ■ 
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